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Quantum fluctuations are of central importance in itinerant ferromagnets; in the case of the Stoner 
Hamiltonian, with contact interactions, they deliver a rich phase diagram featuring a first order 
ferromagnetic transition preempted by a spin spiral and a paired density wave. However, to date 
aU analyses of fluctuation corrections neglect the flnite ranged nature of the Coulomb interaction. 
We develop the formalism to consider the effects of fluctuations with a realistic screened Coulomb 
potential. The finite ranged interaction suppresses the tricritical point temperature of the first order 
ferromagnetic transition, bringing theory into line with experiment, whilst retaining the exotic spin 
spiral and paired density wave. In an ultracold atomic gas a finite ranged interaction damps the 
competing molecular instability, permitting the observation of ferromagnetic correlations. 

PACS numbers: 03.75.Ss, 71.10.Ca, 67.85.-d 



Metallic systems in the vicinity of second order phase 
transitions display remarkable quantum critical phenom- 
ena [T]. However, in several metals [2-6J quantum criti- 
cality gives way to a first order transition, spatially mod- 
ulated magnetic order [THS] , and a p-wave superconduct- 
ing instability [TUHIl]. Starting from a minimal Hamil- 
tonian of freely dispersing electrons with a repulsive con- 
tact interaction, theory indicates that soft magnetic fluc- 
tuations transverse to the magnetic order drive the first 
order transition [T5H2T]. spin spiral phase [15], and the 
p-wave superconducting instability |22| . However, there 
is a significant discrepancy: experiments typically show 
a tricritical temperature of w 0.02rF [51 H HSj (in 
SrgRuaOy), while theory predicts w O.STp [H]. Here 
we demonstrate both analytically and numerically that 
this discrepancy can be resolved by using the realistic 
screened Coulomb potential rather than the ubiquitous 
contact potential. 

A cold atom gas presents an alternative forum to 
explore the Stoner Hamiltonian. Attempts to observe 
itinerant ferromagnetism in a cold atom gas have been 
thwarted piHSO] by a competing instability to a molec- 
ular bound state [31]. Motivated by the suppression of 
pairing in a narrow Feshbach resonance [32J [33] , here we 
demonstrate how a finite interaction range removes pair- 
ing in a region of stable ferromagnetic order. 

In this paper we study itinerant ferromagnetism in the 
presence of a finite ranged repulsive interaction. We 
extend an analytical fluctuation correction formalism 
that has already delivered a phase diagram containing 
a first order ferromagnetic phase transition [16], spiral 
phase [15] , and p-wave superconducting instability |22| 
to now include a finite ranged interaction and calculate 
the phase diagram. Second, we perform complementary 
Quantum Monte Carlo calculations to verify the zero 
temperature behavior. Finite ranged interactions sup- 
press the tricritical point temperature. Finally, we fo- 
cus on the cold atom gas with flnite ranged interactions 



that displays ferromagnetic order whilst circumventing 
the competing pairing process. 



FORMALISM 

To study itinerant ferromagnetism in the presence of a 
finite ranged interaction we focus on the idealized Hamil- 
tonian 

PT.Pi.q 

where — p^/2 — /i is the dispersion, /i is the chemical 
potential, and we adopt atomic units h = m = = \. 
A general momentum dependent interaction g(p,q) acts 
between the two species, where p is the incoming relative 
momentum and q represents the momentum transfer. 
The Stoner Hamiltonian of electrons interacting with a 
contact interaction would be recovered with (/(p, q) = g. 
To study the onset of ferromagnetic order, two ah ini- 
tio, complementary, and mutually consistent tools have 
emerged: the analytical order by disorder approach [16| . 
and Quantum Monte Carlo [TH[33]. We now extend both 
methods to probe and understand the consequences of fi- 
nite range interactions. 



Functional integration 

The functional integral formalism calculates the quan- 
tum partition function expressed as a coherent state field 
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TABLE I. The rescaling of the interaction strength in the 
mean-field and fluctuation contributions to the free energy 
on introducing finite range interactions for the solid state and 
cold atoms cases. The bottom row shows the expected inter- 
action strength of the magnetic transition in the mean-field 
approximation. 
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where the field describes a two component Fermi gas. 
We use the four-momentum notation p = {a;,p} with 
Matsubara frequencies uj and inverse temperature /3 ~ 
l/T. 

We decouple the quartic interaction term with a 
Hubbard-Stratonovich transformation into the full vec- 
tor magnetization and density channel p [16) . This de- 
coupling scheme is sensitive to the spin spiral instability 
driven by transverse magnetic fluctuations, and simul- 
taneously the p-wave superconducting instability driven 
by longitudinal magnetic fluctuations. The action is now 
quadratic in the Fermionic degrees of freedom, and after 
integrating them out we recover the quantum partition 
function Z = J 'D(j)'DpcxTp{—S) with an action 

S = Tr[(l)g(p-pgp]~Tr\n {dr+ia+gp)l- g(p ■ cr , (3) 

where g is the operator form of the potential. We ex- 
pand the action to quadratic order in fluctuations of 
p and around their putative saddle-point values po 
and Mq = M(cos Q • r, sin Q • r, 0), where Q is the spi- 
ral wave vector. A gauge transformation renders the 
magnetization uniform and modifies the dispersion as 
= p2/2 + cr\/(P-Q) + (.9^^)^ - After integrat- 
ing over magnetization and density fluctuations the free 
energy is 

F =X1^ "(Cp<T)+Xl-9(Pt-Pi'0)"(?PTt)"(^Pu) 



Pi. Pi 
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The first term is the kinetic energy, the second the mean- 
field contribution of the interactions, and the third is 
the fluctuation correction. The momentum summation 
in the mean-field term introduces a weighted average 
that can be performed analytically for certain potentials 
at low temperature. The momentum summation in the 
fluctuation term is dominated by the contributions at 
+ = 4fc| [15] and therefore the interaction poten- 
tial can be approximated by ^(p, q) g{^/2kY , ^/^k^) . 
These two observations allow us to simply rescale the 
interaction strength that appears in the mean-field and 
fluctuation correction terms and afterwards treat them 
as pure contact interactions. The validity of the approx- 
imation will be verified in the cold atom section. 

We now enumerate the interaction potentials that we 
adopt for describing the solid state and cold atom gas. 

Solid state: We use the Coulomb interaction 
gQ-r/b I ^^^2^ with screening length 6, whose 
Fourier transform (7(p, q) = 5/(1-1- depends 
only on the momentum transfer q. 

Cold atoms: The T-matrix that describes the Fesh- 
bach resonance can be modeled by several poten- 
tials [m 1311 [35]. Here we concentrate on two: 
the separable form used by Pekker that facili- 
tates the momentum summation [321I33|: 5(p,q) ~ 
{2k-pa / TTv-p)[l + 2ara{p'^ + q^)]] and to establish a di- 
rect contact with the screened Coulomb interaction 
a reciprocal form (a Taylor expansion of the former) 
5(p,q) = (2A:Fa/7rj/F)/[l - 2arc(p^ + g^)]"^ Both 
potentials depend on the energy in the center-of- 
mass frame that appears in the scattering ampli- 
tude. Here the interaction strength is analogous to 
the scattering length k-pa, the effective range is re, 
and z/p is the density of states at the Fermi surface. 

In Table ([l]) we summarize how the interaction potentials 
rescale the mean-field and fluctuation correction poten- 
tials after performing the momentum summations. 



Quantum Monte Carlo 

The fluctuation corrections included in the analytical 
formalism represent a subset of all possible contributions 
to the free energy. To gauge the effectiveness of our ana- 
lytical formalism we perform Diffusion Monte Carlo cal- 
culations with the CASINO program [37]. The approach 
is a refinement of that used in previous studies of itiner- 
ant ferromagnetism [TH[M1[3S]. This method optimizes a 
trial wave function at zero temperature to yield the exact 
ground state energy, subject only to a fixed node approx- 
imation, and thus neatly complements the analytics. We 
use a variational wave function i/; = e^"^ D that is a prod- 
uct of a Slater determinant, D, that takes full account of 
the Fermion statistics, and a Jastrow factor J to include 



3 




0.7 0.8 0.9 1 

9^ 



FIG. 1. (Color online): Phase boundaries at temperatures 
T G {0,0.1, 0.2, 0.305}rF with a screened Coulomb interac- 
tion of strength gv and range fcpfc. The ferromagnetic tran- 
sition is first order (solid blue curves) at low fepfe, separated 
from the second order at larger fcpb (dotted blue curves) by 
the green dots. The dashed red line is the locus of tricriti- 
cal temperatures. The black points with error bars show the 
DMC r = estimates for the phase boundary. The inset 
shows the tricritical temperature (solid black line) and the p- 
wave superconducting critical with fepb, the red dots denote 
the CaRuOa family and ZrZn2. 



further correlations. We used a screened Coulomb repul- 
sion ge~^/^ / Airb'^r that exactly reflects the potential used 
in the analytics. 

The Slater determinant consists of plane-wave spinor 
orbitals containing both spin-up and spin-down electrons, 
D = det({V^kGfeF,t''^kGfcF,4.})- Fixing the spin Fermi sur- 
faces {fep.-f, /cp,^} sets the magnetization. For computa- 
tional efficiency we factorize the Slater determinant into 
an up and a down-spin determinant [37] • Provided that 
the orbitals of the minority spin state are the lowest en- 
ergy orbitals of those in the majority spin state, [39], this 
gives the state whose total spin is S^^^. 

The Jastrow factor, J, accounts for electron-electron 
correlations. It consists of the polynomial and plane- 
wave expansions in electron-electron separation proposed 
in Ref . |40) . To further optimize the wave function the or- 
bitals in the Slater determinant were evaluated at quasi- 
particle positions related to the electrons through a poly- 
nomial backflow function [41] that partially relieves the 
fixed node approximation. The trial wave functions were 
optimized using QMC methods using VMC, backflow, 
and Diffusion Monte Carlo. Twist averaging was em- 
ployed to remove finite size effects. 



PHASE DIAGRAM 

The analytical and computational tools developed cal- 
culate the free energy of the electron gas with screened 
Coulomb interaction. We seek the magnetization that 
minimizes the free energy to construct the phase dia- 
gram. To orient the discussion we first start with the 
common contact interaction dX k^b In Fig. [T] fluctu- 
ation corrections [151 116] drive the transition first order 
at gv w 0.7 (versus the mean-field second order transi- 
tion at gv = 1). The importance of fiuctuations is re- 
duced as we increase the temperature and ultimately the 
transition becomes second order at the tricritical temper- 
ature Tc — O.STp. This tricritical point temperature is in 
agreement with previous studies of the contact interac- 
tion ([151 If6j) but is an order of magnitude higher than 
that seen in experiment. 

At zero temperature, increasing the screening length 
diminishes the interaction strength for the fluctuation 
correction as 1/(1 -I- 4A:|6^) (Table With the driv- 
ing force of the first order transition suppressed, the 
critical interaction strength rises towards the mean-field 
critical interaction strength gv = 1. Nevertheless, the 
transition remains resolutely first order due to a non- 
analytic contribution to the free energy of the form 
M-* log (ff^M^ -h T^) [Tg. However, with mcreasmg k-pb 
a slight temperature rise occludes the non-analyticity and 
the transition reverts to second order, exemplified in the 
inset of Fig. [T] which shows how the tricritical temper- 
ature rapidly reduces with rising screening length. Fo- 
cusing on the CaRuOa family with k^h = 0.40 [H] and 
ZrZn2 with k-ph = 0.47 [43j . the tricritical point temper- 
ature is reduced to O.OSTp and 0.02rF respectively with 
reasonable agreement to experiment [8] H] [23] . 

We also study the magnetic transition at zero tem- 
perature using DMC. This predicts a phase boundary 
with increasing screening length that is quantitatively 
similar to the analytical predictions. Here the Monte 
Carlo is in better agreement with the analytical predic- 
tion compared with previous DMC studies [15l [34l [38] 
that employed a square potential with uncontrolled ef- 
fective screening length. 

A Landau expansion with the zero range system 
demonstrates that the first order ferromagnetic transition 
is always accompanied by a spin spiral phase [13 HI] . To 
incorporate the screening length we take the same Lan- 
dau expansion [22], and rescale the interaction parame- 
ters according to Table (|l]). This reveals that the spiral 
phase persists in a thin strip of order /S.gv « 4 x 10~^ 
between the first order phase boundary and the Lifshitz 
line, terminating at the tricritical point. The transition 
temperature of an instability to a p-wave superconductor 
was calculated by Ref. [35] and the associated interac- 
tion strengths can be rescaled following the prescription 
in Table (|l]). In the inset of Fig. [l] we find that the peak 
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FIG. 2. (Color online) The phase diagram for a T = cold 
atom gas with pseudopotential interaction strength kpa and 
effective range fcpre. The two shaded regions denote zero 
loss (shaded blue), and a loss rate less than O.Sep (shaded 
red). The thick black line represents the transition to the 
ferromagnetic state for the separable pseudopotential (whose 
fluctuation correction rescaling is zero at the black dot), and 
the dotted line the reciprocal pseudopotential transition. The 
inset compares the phase diagram for the full fluctuation cor- 
rections with separable pseudopotential (solid black line) to 
one obtained using the g{p, q) g(\/2fcF, v^fcp) approxima- 
tion (dashed black line). 



superconducting transition drops with screening length 
even more rapidly than the tricritical temperature, ob- 
scuring the instability in typical materials. 



Ultracold atomic gas 

As in the solid state, finite range interactions alter the 
nature of the ferromagnetic transition. At zero effective 
range and zero temperature, the separable and reciprocal 
forms for the pseudopotential are identical so in Fig. [2] 
both magnetic transitions concur with previous studies 
that employed contact interactions [1111^. Sweeping the 
effective range through zero raises the rescaled interac- 
tion strength, reducing the critical interaction strength. 
At the negative effective range kpTe = — l/S/cFa, the scal- 
ing of the fluctuation term for the separable pseudopo- 
tential is zero, recovering the mean-field free energy for 
which the transition occurs at fepo = 57r/7. At more neg- 
ative effective range the fluctuation term is restored, thus 
reducing the critical interaction strength. In the case of 
the reciprocal pseudopotential, at large effective range 
the rescaled interaction strengths scale as l/fcp'^e, lead- 
ing to a suppression of interaction effects and a quench- 
ing of the ferromagnetic transition at fep'^e ^ —0.68. 
Although here the separable and reciprocal pseudopo- 
tentials lead to different phase behavior, this regime is 
not only beyond the effective range approximation where 



higher order terms would be important, but furthermore 
in a regime where we will now determine that losses dom- 
inate. 

In a cold atom gas, to generate the repulsive interaction 
experimentalists exploit a Feshbach resonance between 
the free atoms in the Fermi sea and the bound state of 
two atoms to generate a positive scattering length. How- 
ever, as the bound state is necessarily lower in energy 
there is an instability to pair formation. The loss rate 
has been studied from both the paramagnetic |32) and 
polaronic |33| standpoints. To determine whether the 
ferromagnetic or pairing instability dominates, we follow 
Ref. [32] we start from Eqn. ^ and we decouple this in- 
teraction using a Hubbard-Stratonovich transformation 
into the Feshbach molecule channel Aq = X^k'^k^'^q-kt- 
At leading order the resulting Lagrangian is 



1 , /■d^pnF(Cp-H3^)+nF(Cp-ai)-l^ 



5(0,q) y(27r)3 it^ - Cp+at ^ 



(4) 



c-i(",q) 



After regularizing the above integral [32j , and specializing 
to an attractive potential with finite range To, we obtain 
an expression for the pairing susceptibility C(a;,q) 
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(5) 



The imaginary component of the pairing susceptibility 
pole represents the bound state pairing rate. We find 
generally that the maximal pairing rate occurs for q = 0. 

In Fig. [2] we overlay the magnetic transition with lines 
of equal pairing rate A. Pairing rate reduces with in- 
creasing positive effective range as the molecules become 
more tightly bound. This leads to a window in the phase 
diagram where the system is both magnetized and there 
are no losses. However, at negative effective range the 
losses occur on a time-scale ~ 0.1 ms which is significantly 
shorter than the trap crossing time ^ 1 ms, so large mag- 
netic domains cannot be formed. To date all Fermionic 
mixtures used in cold atom gas experiments have neg- 
ative effective ranges so are not suitable for observing 
magnetic correlations |44j . however a polar molecule gas 
with strong dipolar interactions does display large posi- 
tive effective range [45] and a positive s-wave scattering 
length and so presents an opportunity to observe ferro- 
magnetic phenomena. 

The inset of Fig. [2] confirms that the phase boundary 
obtained by rescaling the interaction strength of the fluc- 
tuation contribution with the approximation \-\-2arf.{p^+ 
q^) n- 1 + Skpakpr^ conforms with the phase boundary 
resulting from the exact momentum summation. This 
verifies the rescaling approximations given in Table (m|. 
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DISCUSSION 

The analytical fluctuation correction formalism with 
a contact repulsion has successfully demonstrated how 
quantum fluctuations drive not only a first order ferro- 
magnetic transition, but also a spin spiral and a p-wave 
superconducting instability. The finite ranged interac- 
tions reduce the fiuctuation corrections so dramatically 
suppress the tricritical point temperature from O.STp to 
0.02rF. This is quantitatively similar to the w 0.02Tf 
seen in ZrZn2 [SI El [23] . Finite ranged interactions sup- 
press the transverse magnetic fluctuations that occluded 
quantum criticality. With the tricritical temperature 
now reduced this system allows for the observation of 
quantum critical phenomena, such as those envisaged by 
Hertz. 

The idealized Stoner Hamiltonian can also be studied 
in a cold atom gas. The positive effective finite range 
interaction, found in a polar molecular gas, eliminates 
the pairing permitting the observation of ferromagnetic 
order. 
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